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Abstract 

Data from electron-proton collisions at a center-of-mass energy of 300 GeV are used 
for a search for selectrons and squarks within the framework of the minimal super- 
symmetric model. The decays of selectrons and squarks into the lightest supersym- 
metric particle lead to final states with an electron and hadrons accompanied by 
large missing energy and transverse momentum. No signal is found and new bounds 
on the existence of these particles are derived. At 95% confidence level the excluded 
region extends to 65 GeV for selectron and squark masses, and to 40 GeV for the 
mass of the lightest supersymmetric particle. 
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1 Introduction 



Supersymmetry is presently considered to be a promising candidate for a theory beyond the 
Standard Model (SM) of particle physics. In particular the Minimal Supersymmetric extension 
of the Standard Model (MSSM) describes as well as the SM all experimental data, and in 
addition it offers solutions for some of the questions left open by the SM, such as the Higgs mass 
or the hierarchy problem. However, no direct evidence for supersymmetry has yet been found. 

Supersymmetry relates fermions to bosons and predicts for each SM particle a partner with 
spin differing by half a unit. So sneutrinos and selectrons, {i'eL,eL), sr, are scalar partners of 
neutrinos and electrons, (i^eL, e^), e/?, and similarly squarks {uL,dL),uji,dji are the partners of 
up and down quarks. Two Higgs doublets with vacuum expectation values f2,fi are necessary 
to give masses to up-type quarks (^2) and to down- type quarks and charged leptons (vi). The 
partners of the gauge bosons W^,Z^,-) and the two Higgs doublets are called gauginos and 
higgsinos. They can mix and form two charged mass eigenstates X12 (charginos) and four 
neutral mass eigenstates X1234 (neutralinos). 

Experimentally supersymmetric particles are constrained to be heavier than their SM part- 
ners and so supersymmetry must be broken. In the MSSM this leads to extra mass parameters 
M2 and Ml for the SU2 and Ui gauginos. Thus the masses of charginos and neutralinos depend 
on Ml, M2, tan/? = v^jvx and the higgsino mass parameter [i. The gaugino-sfermion couplings 
are the same as in the SM, but via mixing the chargino and neutralino-sfermion couplings also 
depend on these supersymmetric parameters. 

In the MSSM it is assumed that the multiplicative quantum number i?-parity is conserved, 
where Rp = 1 for the SM particles and Rp = —1 for their superpartners. This implies that 
supersymmetric particles can only be produced in pairs and that the lightest of them, which 
generally is assumed to be the Xi; is stable. Since the Xi is neutral and only weakly interacting 
it will escape direct experimental detection. 

At HERA the dominant MSSM process is the production of a selectron and a squark via 
neutralino exchange ep — > e g X as shown in Fig. The e and q decay into any lighter gaugino 
and their SM partners. The decay involving gives an experimentally clean signature of an 
electron, hadrons and missing energy and momentum. This is the channel analyzed here, where 
the proper branching ratios for e — > ex? and q qxi are taken into account. 




Figure 1: Feynman diagram for selectron - squark production via neutralino exchange and the 
subsequent decays into the lightest supersymmetric particle Xi- 

In this paper it is assumed that the masses of and ep_ are identical. The same holds 
for the q states and mass degeneracy is also assumed for the superpartners of the four lightest 



quarks. In order to reduce the number of free MSSM parameters it is assumed that Mi and 
M2 are related via the weak mixing ang le 9w, Ml = 5/3 M2 tan^^v^, as suggested by Grand 
Unified Theories (GUT's) No other GUT relations are used. The parameter set given above 
specifies completely all couplings and masses involved. 

At e~^e~ experiments masses below 45 GeV for selectrons and for squarks have been excluded 
by direct searches and no attempt was made here to evaluate this excluded region. A search 
for selectron production at 130 and 136 GeV ^ yielded a lower limit of Mg = 56 GeV for 
M^o < 35 GeV, but however, with different assumptions for the ept, bl mass splitting. 

At the resonance direct searches for charginos and neutralinos H] restrict the MSSM 
parameter space mainly for higgsino-like Xi ^or Xi masses up to 20 GeV. These results 
are independent of Mg and Mi>. Above the resonance recent searches have extended these 
limits considerably ^. Here also e and 9 exchange contributes and, due to destructive 

interference, the constraints are most stringent for large Mg and Mj^. 

The experiments at the Tevatron pp collider have obtained strong bounds on q masses. These 
limits however either depend on the assumption of a light Xi or they are only valid for small 
gluino masses [^. The latter bounds can only be related to the results obtained here if one 
assumes additional GUT relations between the gluino mass and M2, Mi. 

While this search deals with supersymmetric models with conserved Rp, the production of 
q in models with broken Rp has been considered in Q. 



2 Detector description 



A detailed description of the HI apparatus can be found elsewhere [10|. The following briefly 
describes the components of the detector relevant to this analysis, which makes use mainly of 
the calorimeters. 



The liquid Argon (LAr) calorimeter ||ll| extends with full azimuthal coverage over the polar 
angular range 4° < ^ < 153°, where 6 is defined with respect to the proton beam direction (+z 
axis). The calorimeter is highly segmented and consists of an electromagnetic section with lead 
absorbers, corresponding to a depth of between 20 and 30 radiation lengths, and a hadronic 
section with steel absorbers. The total depth of the LAr calorimeter varies between 4.5 and 8 
hadronic interaction lengths. Single particles are measured in the LAr calorimeter with energy 
resolutions ofa{E)/E ^ 0.11/v/^/GeVeO.Ol for electrons and a{E)/E 0.5/y^/GeVe0.02 
for charged pions. The electromagnetic and hadronic energy scales are known to 2% and 5% 
respectively. 

The calorimeter is surrounded by a superconducting solenoid providing a uniform magnetic 
field of 1.15 T parallel to the beam axis in the tracking region. A large instrumented iron 
structure surrounds the solenoid to serve as a return yoke for the magnetic flux and as an 
additional calorimeter of 4.5 hadronic interaction lengths to absorb tails of hadronic showers. 

The backward region of the detector (151° < < 177°) was covered in 1994 by an electro- 
magnetic lead-scintillator calorimeter, which was replaced in 1995 by a new lead-scintillating 
fibre calorimeter [^] with improved containment and angular coverage. A plug calorimeter 
covers the angular region 0.75° < 6 < 3.5° around the forward beam pipe. 

Charged particle tracks are measured in two concentric jet drift chambers and a forward 
tracking detector covering the polar angular range 7° < 6 < 165°. In this analysis the tracking 
detectors are used for the determination of the event vertex only. 



3 Data selection and analysis 



The event sample corresponds to an integrated luminosity of 6.38 it 0.13 pb~ and was collected 
by the HI experiment in 1994 and 1995 during collisions of 27.5 GeV positrons^ with 820 GeV 
protons. 

First some event quantities which are used in the analysis are introduced. The quantity E—P^ 
is defined as -E'i(l — cos Oi) where the sum runs over all LAr and instrumented iron calorimeter 
energy deposits. The backward and plug calorimeters are used only as a veto. The (2 dimen- 
sional) transverse momentum vector of a single energy deposition is Pt^i = Ei sin 0j(cos sin (pi) 
with polar angle 9i and azimuthal angle 4>i. The total transverse momentum vector Pt is the 
vector sum of all Pt^i in the LAr and iron calorimeter. The variable Pt\\ {Pt±) is defined as 
the absolute value of the component of Pt which is parallel (perpendicular) to the transverse 
component of the identified scattered electron. 

The following selection criteria are used: 

1. Events must be accepted by a LAr calorimeter trigger for electrons which requires a local 
energy deposit of more than 8 GeV in a small region of the electromagnetic LAr calorimeter. 

2. A vertex has to be found from measured charged tracks to be within 35 cm of the nominal 
interaction point. 

3. An electron must be identified by its shower shape in the LAr calorimeter, it must have 
an energy greater than 10 GeV, and its transverse momentum has to exceed 8 GeV. The 
polar angle 9e must be in the range 10° < 0^ < 135° and the electron must be isolated in 
a cone of pseudorapidity and azimuthal angle of radius 0.5. 

4. The energy deposited in the backward calorimeter has to be below 5 GeV. 

5. The transverse momentum measured in the plug calorimeter must be below 2 GeV. 

6. Filters against cosmic and beam halo muons are applied. 

7. The quantity E — P^ must be below 40 GeV. 

8. Either 

(a) the polar angle of the "current quark" must fulfill 9q > 5°. Assuming standard 
deep inelastic scattering kinematics, Oq can be calculated in the quark parton model 
from the scattered electron and the electron beam energy. To compensate for an 
undetected radiative photon in the initial state, the electron beam energy is reduced 
by the energy of this photon which is determined using the measured total E — Pz- 

(b) Pt± > 3 GeV and 

(c) Pt\\ > 2.5 GeV+ - Pz)/2 

or 

(d) Pt± > 7 GeV and 

(e) > 3 GeV 

The asymmetry of the HERA beam energies leads to a boost in the proton beam direction of 
the produced heavy e and q and their decay products. Therefore, in contrast to the majority 
of SM events, no energy deposition is expected in the backward region of the detector (cuts 3 
and 4). Energy losses in the forward beam pipe can create large values of Pt- Such events are 
vetoed by cut 5. 

^In this paper the generic name "electron" is used also for positrons. 



Because of momentum conservation, the value E — Pz = 2 x 27.5 GeV = 55 GeV of the 
initial state beam particles is expected also for the measured final state of most SM processes. 
Supersymmetric events would have a much smaller value of E — due to the undetected Xi's- 
This is demonstrated by the dashed histogram in Fig. ||a which shows a superposition of several 
hundred simulated e q events with Mg 



65 GeV and various values of M^o. 
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Figure 2: Comparison of data (points) with the expectation for all contributing SM processes 
(solid histogram). To indicate the signal region, simulated e q events are shown with arbitrary 
normalization by the dashed histogram for Mg = Mq = 65 GeV and M^o between 35 and 55 
GeV. Figure a) shows the E — P^ distribution after basic electron selection cuts 1-6. Figures b) 
and c) show the Pt component perpendicular (Ptj_) and parallel (Pt\\) to the electron direction 
after applying an additional cut of E — P^ < 40 GeV. Figure d) shows data (solid dots) and 
simulated e q events (open circles) after applying cuts 1-7, 8a and 8b (see text). The combined 
cut on > 2.5 GeV + {E — Pz)/2 is indicated by the dashed line. 



The event generator HERASUSY [13|, which is used for this simulation, is based on the cross 
section for all four neutralino exchange diagrams and their interference |14]. All couplings and 
masses involved are determined using a modified version of the ISASUSY program [^]. The 



PYTHIA program [^] is used for parton showers, string fragmentation and particle decays and 
a detailed simulation of the HI detector response is included. 

The events selected by cuts 1-6 are shown as solid points with statistical errors in Fig. 
They exhibit a clear peak at E — = 55 GeV and a long tail towards lower values. This 
distribution is well described by the solid histogram which represents a simulation of the con- 
tributing SM processes, i.e. neutral current and charged current deep inelastic scattering events 
simulated using the DJANGO program and photoproduction events simulated using the 
PYTHIA program |]l6|. The tail towards low E — is due to particles lost in the —z direction 
and is populated mainly by neutral current events with a photon radiated off the initial state 
electron. Contributing at very low E — P^ are also photoproduction processes, in which the 
scattered electron is not detected and a hard photon is misidentified as an electron in the event 
selection. 

After cut 7 motivated by the HERASUSY simulation Fig. ^b and c show the distributions 
of Pt± and P^ii for the remaining 384 events. 367 it 19 events are expected from all background 
sources. The resolution for Pt± is much better than for P^y since more energy is deposited 
collinear to the electron in SM processes. The extreme tails of the Pt distributions are populated 
by radiative charged current events with a photon misidentified as an electron. In all three 
distributions no obvious excess over expectation is observed. 

Cut 8a on the jet angle predicted by the measured electron is employed in order to reduce 
the neutral current background of events with jets lost in the beam pipe. Also a minimum value 
of Pt± is required (cut 8b). Cut 8c on the correlation of Pj|| and E — P^ is illustrated in Fig. Qd 
which shows the data after cuts 1-7 and 8a together with a sample of simulated e q events. Cut 
8c leads to a considerable loss for e q events with small values of . A large fraction of these 
events are recovered by an alternative choice of cuts requesting mainly a large value of Pt± (cuts 
8d and 8e). 

After all cuts no events are found in a simulation of neutral current (photoproduction) 
processes in which there were 3 (0.63) times the statistics of the data. A background of 0.6 it 0.2 
radiative charged current events is expected. No candidate event remains in the data. 



4 Results 

Having observed no signal, exclusion limits for e q production are derived. 

The selection efficiency is determined using events generated as described above with the 
HERASUSY program for different values of Mg, Mq and M^o. Other MSSM parameters do not 
influence the final state kinematics significantly. In order to interpolate between the simulated 
values for Mg, Mq and M^o an empirical function e('P) is fitted to the simulated efficiency values. 
Here V is defined via 

2 Me 2 Mq 

It characterizes the transverse momenta of electron and quark after the 2-body decays e — s- exi 
and q —* qxi- The efficiency €{V) reaches a plateau of 56% for V > 25 GeV. For smaller values 
of V any signal in the detector becomes less and less significant and the efficiency drops to 50%, 
18% and 0% for values of P = 20, 10 and 5 GeV respectively. 

Systematic errors of this analysis originate in the efficiency parameterization, the electromag- 
netic and hadronic calorimeter energy scales, the parton density functions and the integrated 



luminosity. The total systematic uncertainty amounts to 4.3% in the number of expected su- 
persymmetric events. 



In order to determine excluded regions in the MSSM parameter space, a systematic scan 
is performed for values of Mg > 45 GeV, Mg > 45 GeV, -1 TeV < ^ < 1 TeV, Ma > GeV 
and 1 < tan/3 < 50. For each chosen set of these parameters the cross section, branching ratios 
and efficiency are evaluated yielding an expectation for the number of observed signal events. A 
parameter set is excluded when this number exceeds the 95% confidence level upper limit for 
observed events taking into account the 4.3% systematic uncertainty mentioned above. Details 



of this statistical procedure are described in |18]. The resulting excluded regions in the MSSM 
parameter space are presented in the planes 

• M^o versus (Mg + Mq)/2 for constant fi and tan/3 = 1.41 (Fig. |) 

• Me versus Mg for constant M2, ^ and tan/3 = 1.41 (Fig. ^ 

• M2 versus /x for constant (Mg + Mg)/2 and tan/3 = 1.41 (Fig. H) 

• M2 versus tan/3 for constant (Mg + Mq)/2 and /i (Fig. |^) 



The choice of tan/3 = 1.41 is taken as typical of a low value for tan/3, and allows direct comparison 
with recent results obtained at LEP. The dependence on tan/3 is shown in Fig. ^. 

The cross section for e q production depends to first approximation on the Xi mass and 
couplings, and on the energy threshold for the hard process Mg + Mq. The neutralino and 
chargino masses scale roughly with M2 for /i ^ 0. For the excluded parameter regions the 
branching ratios B(e Xi) ^{q Xi) ^'^^ close to one because all other gauginos are 
heavier than Afg or Mq. Only at very small jj, and M2 other decays are possible, which have not 
been searched for in this paper. 

In Fig. m the plane M^o versus (Mg + Mg)/2 is shown for tan/3 = 1.41. The dark and light 
shaded areas are excluded for /i = —50 GeV. For smaller values of (^ ^ — cxd) and for small 
M2 the competing decays into other gauginos lead to reduced limits for (Mg + Mq)/2. The 
dark shaded area is excluded for all /i < —50 GeV. In the range /i < —50 GeV, and for this 
fixed value of tan/3, the Xi mass depends on M2 only, which is indicated by the second vertical 
scale on the right side of the figure. At large M^o the mass difference between e, q and y? 
becomes small and the efficiency drops correspondingly. Also the cross section is reduced due 
to the higher propagator mass. Both effects lead to reduced limits for (Mg + Mg)/2. As seen 
from the distance between the shaded area and the diagonal line for M^o = Mg^g this search is 
sensitive down to mass differences of ~ 10 GeV. The excluded mass range extend to 65 GeV 
for (Mg + Mg)/2 and to 40 GeV for M^o. In particular the region (Mg + M^)/2 < 63 GeV is 
excluded for M,,o < 35 GeV. Also shown are the limits obtained from direct searches for e and 
q pair production in decays 

Fig. Q shows the excluded regions in Mg versus Mq for three different values of M2. For a 
large range of Mg and Mq the limits are approximately a function of only the sum Mg + Mg, as 
used in Fig. |3[ This approximation is valid for large mass differences between e and Xi and also 
between q and Xii i-e- Me,q — M^o > 20 GeV. The limits on Mg and Mq decrease when these 
differences become too small. For M2 ~ 50 GeV the highest mass bounds are reached, namely 
Mg = 77 GeV at Mq = 50 GeV and Mq = 75 GeV at Mg = 50 GeV. 

In Fig. ^ the exclusion limits are shown in the M2 versus /i plane for tan/3 = 1.41 and the mass 
ranges (Mg + Mg)/2 < 60 GeV and < 63 GeV. For ^ <C the x? is dominated by its photino 
component. Consequently the coupling of e and q to the Xi are of electromagnetic strength and 
allow for a sizeable cross section. When fi approaches the Xi becomes higgsino-like and the 
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Figure 3: Excluded regions at 95% confidence level in the plane of M^o versus half the sum of 
Mg and Mg for tan/3 = 1.41. The dark and light shaded areas are excluded for /i = —50 GeV. 
The dark shaded area is excluded for all /x < — 50 GeV. The difference between both is due to 
decays into other gauginos than the Xi- The diagonal line corresponds to Mg^g = M^o. The 
hatched area is excluded by direct searches for e and q pair production in decays The 
vertical axis is also labeled in terms 0/M2 (see text). 
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Figure 4: Exclusion limits at 95% confidence level in the Mg versus Mg plane for different 
values of M2 and /i < —50 GeV, tan/3 = 1.41. Mg and Mg masses below 45 GeV are excluded 
by searches in decays. 
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Figure 5: Domains excluded at 95% confidence level in the M2 versus /x plane from the search 
fore q production for tan(3 = 1.41. The dark shaded area corresponds to {Me+Mg)/2 < 63 GeV, 
the additional light shaded area to {Mg + Mq)/2 < 60 GeV. The hatched areas show the most 
restrictive bounds obtained from e~^e~ collisions at LEP at center- of -mass energies around the 
mass and at 130 and 136 GeV /Q/, see text. 



couplings and cross sections become very small. This defines the limit of the excluded parameter 
space for /x close to 0. For <C the range excluded for M2 becomes independent of /x. The 
region of /i ^ and small M2 is not excluded because of decays into charginos and corresponds 
to the light shaded area in Fig. ^. 

Also displayed as hatched areas in Fig. ^ are the most stringent bounds from searches for 
charginos and neutralinos at LEP at center-of-mass energies around the mass and at 130 and 
136 GeV 1^, ^. These limits are valid for M^^u > 500 GeV. For smaller masses the bounds are 
less restrictive. For ^ <C the search presented here extends to considerably larger values for 
M2 than the LEP results. For ^ > however the LEP results are more restrictive. Therefore 
only values for /i < are displayed. 

In Fig. |6| the dependence of the M2 limits on tan/3 is displayed for (Mg + Mg) jl < 60 GeV 
and = —200 GeV. The limit on M2 rises to 80 GeV for large values of tan/?. The two symbols 
at tan/3 = 1.41 and 35 indicated the bounds from LEP as discussed above p. The range of M2 
excluded in addition to those of LEP is largest for small values of tan/?. 

5 Conclusion 

Within the MSSM with conserved i?-parity we have searched for e q pairs and their dominant 
decay modes into the lightest supersymmetric particle Xi- No signal was found and rejection 




Figure 6: Excluded region at 95% confidence level for M2 as a function of tan/3 for 
(Me + Mq)/2 < 60 GeV and fj, = —200 GeV (shaded area). The two symbols show the upper 
limits on M2 obtained for tan/3 = 1.41 and 35 from e^e~ collisions at LEP for center- of -mass 
energies of 130 and 136 GeV /§/. 

limits at 95% confidence level are derived for the parameters Mg, Mg, M2, ^ and tan/?. The 
cross section, the branching ratios and the efficiency depend mainly on the sum Mg + Mq and 
on the Xi mass and couplings to fermions. The search is therefore most sensitive for a x? with 
a large photino or Z component. For tan/3 = 1.41 masses up to (Mg + Mq)/2 = 65 GeV and 
up to M^o = 40 GeV are probed. The dependence on tan/3 is weak and mass limits improve 
towards large tan/3. The parameter region 

(Mg + Mg)/2 < 63 GeV for M^o < 35 GeV and /i = -50 GeV 

is excluded. This limit applies as long as Mg g — M o > 20 GeV and decreases slightly when this 
difference becomes smaller. The result becomes independent of ^ for < —50 GeV. At small 
values of M2 the sensitivity is reduced due to decays into charginos. 

For a photino-like Xi and small values of tan/? these results extend considerably the limits 
obtained by the searches for supersymmetry at LEP 1 and also those obtained while LEP was 
running at center of mass energies of 130 and 136 GeV. The limits, in particular for squarks, 
are independent of the gluino mass. 
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